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Generalized  algorithm  to  build  descriptor  matrices  of  convex  politope  Trefftz  finite  elements  is 
introduced  in  the  paper.  A finite  number  of  plane  waves  is  used  to  expand  electromagnetic  field 
inside  the  elements.  Projections  of  the  intra-element  field  on  an  additional  set  of  basis  functions 
defined  on  the  element  surface  are  used  to  build  admittance  matrix  descriptors  of  the  element. 

INTRODUCTION 

As  a generalization  of  the  method  of  minimum  autonomous  blocks,  introduced  by  V.V. 
Nikol’skii  and  T.I.  Lavrova  in  the  late  70s  [1],  Trefftz  Finite  Element  method  (TFEM) 
has  been  recently  introduced  into  the  computational  electromagnetics  [2],  [3],  As  the 
conventional  finite  element  method  (FEM),  the  TFEM  is  based  on  a division  or 
decomposition  of  a boundary  value  problem  for  Maxwell's  equations  into  a set  of 
elements.  Though,  instead  of  polynomial  functions,  plane-wave  solutions  of  the 
Maxwell’s  equations  are  used  as  the  intra-element  basis  functions  to  expand  electric  and 
magnetic  fields  inside  the  elements.  Though,  the  basic  concepts  of  the  method  are  quite 
general  and  are  outlined  in  [2],  [3],  there  is  no  formalized  procedure  to  build  descriptors 
of  complex  polytope  structures  such  as  polygonal  prisms,  tetrahedrons  and  so  on.  This 
paper  introduces  such  formalized  procedure  to  build  the  admittance  matrices  of  convex 
polytope  elements. 


BUILDING  DESCRIPTORS  OF  3-D  ELEMENTS 


Let  us  consider  a 3-D  boundary  value  electromagnetic  problem  in  the  frequency 
domain.  The  problem  is  described  by  the  Maxwell’s  equations  and  boundary  conditions 
in  a Cartesian  coordinate  system.  The  problem  is  subdivided  into  a set  of  small  convex 
polytope  elements.  All  external  and  internal  boundaries  of  the  problem  are  mapped  on 
the  boundaries  of  the  elements.  An  element  can  be  represented  as  a convex  polytope 
Qp  in  three-dimensional  Euclidian  space  with  N fMr  polygonal  faces  Fn,n  = 1,...,  N fM,t,  . 

The  element  is  uniformly  filled  with  an  isotropic  medium.  Let  us  expand  the  polytope 
element  interior  field  using  Nml  pairs  of  plane  wave  solutions  of  the  Maxwell’s 


equations.  The  field  distribution  inside  the  element  can  be  expressed  as 


'm' 

H{r\ 


= 1 


a: 


F 

77 


+ A.., 


(i  in  j 


(E  \ 

^0  m 1 
V " 0m  J 


ren,k*E<lm=0JI„  = 


7T  k.„  X E„ 


Zn 


(1) 


where  symbol  * denotes  scalar  products,  symbol  x denotes  vector  product,  km  is  the 
unit  vector  of  propagation  direction  of  the  plane  wave  number  m,  r is  a radius  vector  of 
a point  inside  the  polytope  or  on  its  boundary,  E0m  is  the  unit  electric  field  vector  of  the 


plane  wave,  H0m  is  the  magnetic  field  vector  of  the  wave,  and  Am  are  unknown 


magnitudes  of  the  waves  or  expansion  coefficients,  k0  is  a propagation  constant  of  the 


plane  wave,  and  Z0  is  a characteristic  impedance  of  the  plane  wave: 
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k = co-JUju,  Z0  = yjfi/e  . (2) 

Nint  magnitudes  of  the  plane  waves  in  (1)  are  considered  to  be  linearly  independent. 

To  provide  connectivity  of  the  element  interior  with  the  surrounding  elements  and  to 
impose  the  inter-element  boundary  conditions,  an  additional  set  of  basis  functions  is 
defined  on  the  faces  of  the  polytope  element.  Let  us  choose  Ln  vector  basis  functions 

defined  on  the  polygonal  face  Fn  for  the  electric  field  expansion  and  Ln  basis  functions 
on  the  same  face  to  expand  the  magnetic  field.  The  total  number  of  pairs  of  the  surface 
field  expansion  functions  is 

N face 

N,„  = Y.L,  0) 

n= 1 

Designating  the  surface  basis  functions  defined  on  the  face  Fn  as  en(l)  and  hn(l) , the 

electric  and  magnetic  fields  on  the  surface  of  the  polytope  element  can  be  expressed  as 
follows: 

_ I;,  _ Njace  L„  _ 

^ surf  ~ Vn(l)  ' en(l)  > H surf  = XI /L  W)  ’ ^«(0  (^) 

n=l  /= 1 n= 1 /=1 

where  vn(/),  iM(/)  are  unknown  expansion  coefficients.  The  total  number  of  the  boundary 


field  expansion  coefficients  in  (4)  is  2Nsurf . It  is  in  addition  to  the  Nmx  independent 


interior  expansion  coefficients  (1).  The  number  of  the  interior  basis  functions  Nmt  can 
be  chosen  equal  to  the  number  of  the  electric  or  magnetic  field  surface  basis  functions 
Nsurf . It  provides  a possibility  to  uniquely  define  a matrix  descriptor  of  the  element.  To 

do  so,  we  can  project  the  interior  field  on  the  surface  basis  functions  [2],  [3]: 

Vn(l)  K(l)  ~ ^n(l)  }’  / — 1 17  — I,...,  N j-ace  (5) 

where  E and//  are  values  of  the  electric  and  magnetic  fields  defined  by  (1)  taken  on  the 
face  Fn  where  basis  functions  en(l)  and  hll{l)  are  defined.  Either  point  matching  or 

Galerkin  projectors  can  be  used  in  (5).  Constant  vector  basis  functions  may  be  used  with 
the  point  matching  projectors  defined  as 

p«{A^(/)}=D(fn(/))*6„(/)(i;,(/)),  / = 1,...,T„,  n = face  (6) 

where  D is  either  E or  // , bn(l)  is  either  en{l)  or  hn(l),  rn(l)  is  the  radius  vector  of  a 


matching  point  on  the  face  Fn  for  the  basis  function  /.  The  matching  points  can  be 

defined  as  centroids  of  the  polygonal  areas  where  the  corresponding  face  basis  functions 
are  defined.  Galerkin  or  averaging  projectors  can  be  defined  as 

p« {a^(,)}=  — !— | \D*b*(,yds,  l = 1 ,...,Ln,  n = 1 ,...,Nface  (7) 


where  the  integral  is  taken  over  the  surface  of  Fn , and 


Nn{l)  are  the  norms  of  the 


expansion  functions.  Substituting  (1)  into  (5)  we  can  obtain  the  following  relations 
between  the  interior  field  expansion  coefficients  and  the  surface  expansion  coefficients: 
v=M^-A++M;-A~,  J = m; -A+  +M~h  -A-,  (8) 
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where  A + and  A ‘ are  vectors  with  Nmt  coefficients  A * and  Am  ( 1 ),  / and  v are 


vectors  with  Nsurf  components  defined  as 

1 = fe’-- v, j-  f ’ V = J-  A = [v„(1),...,v, , (9) 

where  symbol ' denotes  transposition. 

Matrices  Mc:  and  Mft  are  N surf  by  Nml  complex  matrices  of  projections  of  the  interior 

basis  functions  (1)  on  the  boundary  basis  functions  with  the  elements  defined  with 
either  projectors  (6)  or  (7)  as 
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m = I = \,...L,rn  = \,...,Nfaci, 

A descriptor  matrix  of  the  polygonal  element  can  be  deduced  by  eliminating  unknown 
interior  field  expansion  coefficients  A 1 from  (8).  From  here  on  we  assume  that 
N im  = Nxlirf  = N , which  leads  to  square  N by  TV  matrices  andA/^  . An  assumption 


of  the  equality  of  two  terms  in  the  magnetic  field  projection  sum  or  alternatively  in  the 
electric  field  projection  sum  (8)  leads  to  two  alternative  additional  expressions  to 
construct  the  descriptor  matrix: 

A-=T,-A\  Th=(M-y'-M ; (11) 

A"  = 7j,  • A + , Te=(M;y'-M;  (12) 

Now,  an  admittance  matrix  descriptor  of  the  element  relating  the  unknown  boundary 
expansion  coefficients  can  be  defined  as 

J = Y -v,  Y = Mh  ■ M~',  YeC AVV,  (13) 

where  Me  = M;  + M]  • Tclh,  Mh  = M;  + Mh  • T,:h. 

The  linear  independency  of  the  element  interior  basis  functions  is  the  necessary 
condition  of  existence  of  a non-degenerate  descriptor  of  a polytope  element.  Plane 
waves  propagating  in  the  directions  perpendicular  to  the  sides  of  a convex  polytope 
provide  such  a system  of  functions.  Assembling  of  the  admittance  matrices  (13)  into  a 
global  admittance  matrix  is  a simple  and  straightforward  procedure  and  is  described  in 
[2j. 

CONCLUSION 


Trefftz  finite  element  method  has  been  generalized  in  the  paper  on  the  problems 
subdivided  into  a set  of  convex  polytope  elements.  The  boundary  value  problem  has 
been  reduced  to  a building  and  re-composition  of  admittance  matrices  of  the  polytope 
elements.  Generalized  matrix  formulas  are  derived  to  build  the  admittance  matrix 
descriptors  of  the  elements. 
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